On Optimal Finite-Difference Approximation of PML
نویسندگان
چکیده
A technique derived from two related methods suggested earlier by some of the authors for optimization of finite-difference grids and absorbing boundary conditions is applied to discretization of perfectly matched layer (PML) absorbing boundary conditions for wave equations in Cartesian coordinates. We formulate simple sufficient conditions for optimality and implement them. It is found that the minimal error can be achieved using pure imaginary coordinate stretching. As such, the PML discretization is algebraically equivalent to the rational approximation of the square root on [0, 1] conventionally used for approximate absorbing boundary conditions. We present optimal solutions for two cost functions, with exponential (and exponential of the square root) rates of convergence with respect to the number of the discrete PML layers using a second order finitedifference scheme with optimal grids. Results of numerical calculations are presented.
منابع مشابه
An optimal 25-point finite difference scheme for the Helmholtz equation with PML
In this paper, we analyze the defect of the rotated 9-point finite difference scheme, and present an optimal 9-point finite difference scheme for the Helmholtz equation with perfectly matched layer (PML) in two dimensional domain. For this method, we give an error analysis for the numerical wavenumber’s approximation of the exact wavenumber. Moreover, based on minimizing the numerical dispersio...
متن کاملOptimal order finite element approximation for a hyperbolic integro-differential equation
Semidiscrete finite element approximation of a hyperbolic type integro-differential equation is studied. The model problem is treated as the wave equation which is perturbed with a memory term. Stability estimates are obtained for a slightly more general problem. These, based on energy method, are used to prove optimal order a priori error estimates.
متن کاملAnalysis of a finite element PML approximation for the three dimensional time-harmonic Maxwell problem
In our paper [Math. Comp. 76, 2007, 597–614] we considered the acoustic and electromagnetic scattering problems in three spatial dimensions. In particular, we studied a perfectly matched layer (PML) approximation to an electromagnetic scattering problem. We demonstrated both the solvability of the continuous PML approximations and the exponential convergence of the resulting solution to the sol...
متن کاملAPPROXIMATION OF STOCHASTIC PARABOLIC DIFFERENTIAL EQUATIONS WITH TWO DIFFERENT FINITE DIFFERENCE SCHEMES
We focus on the use of two stable and accurate explicit finite difference schemes in order to approximate the solution of stochastic partial differential equations of It¨o type, in particular, parabolic equations. The main properties of these deterministic difference methods, i.e., convergence, consistency, and stability, are separately developed for the stochastic cases.
متن کاملAnalysis of a finite PML approximation to the three dimensional elastic wave scattering problem
We consider the application of a perfectly matched layer (PML) technique to approximate solutions to the elastic wave scattering problem in the frequency domain. The PML is viewed as a complex coordinate shift in spherical coordinates which leads to a variable complex coefficient equation for the displacement vector posed on an infinite domain (the complement of the scatterer). The rapid decay ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Numerical Analysis
دوره 41 شماره
صفحات -
تاریخ انتشار 2003